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ABSTRACT 

We examine the effects on the Fundamental Plane (FP) of structural departures 
from an galaxy light profile. We also explore the use of spatial (i.e. volumetric) 
as well as projected galaxy parameters. We fit the Sersic law to the V-band 
light profiles of 26 E/SO Virgo galaxies, where n is a shape parameter that allows for 
structural differences amongst the profiles. The galaxy light profiles show a trend of 
systematic departures from a de Vaucouleurs law, in the sense that n increases 
with increasing effective half-light radius Re- This results in R^, and the associated 
mean surface brightness within this radius, having systematic biases when constructed 
using an R^^^ law. Adjustments to the measured velocity dispersion are also made, 
based upon the theoretical velocity dispersion profile shapes of the different light 
profiles, constructed assuming spherical symmetry and isotropic pressure support. 

We construct the FP for the case when structural homology is assumed (specif¬ 
ically, an law is fitted to all galaxies) and central velocity dispersions, tJoi are 
used. The plane we obtain is Re oc (ji'33±o.iOj^-o.79±o.ii^ where Ee is the mean surface 
brightness within the projected effective radius Re- This agrees with the FP obtained 
by others and departs from the virial theorem expectation R oc We find that 

allowing for broken structural homology through fitting profiles (with n a free 

parameter), but still using central velocity dispersions, actually increases the depar¬ 
ture of the observed FP from the virial plane — the increase in effective radius with 
galaxy luminosity (and n) is over-balanced by an associated decrease in the mean 
surface brightness. 

In examining the use of spatial quantities and allowing for the different velocity 
dispersion profiles corresponding to the observed light profiles, we find that use of the 
spatial velocity dispersion at the spatial half-light radius decreased the departure of 
the observed FP from the virial plane. (Use of the spatial half-light radius and mean 
surface brightness term had no effect on the FP as they are constant multiples of their 
projected values). Through use of the Jeans hydrodynamical equation, we convert the 
projected central aperture velocity dispersion, ao, into the infinite aperture velocity 
dispersion, tJtot.n (which is equal to one-third of the virial velocity dispersion). Using 
both the fit and crtot.n we obtain i?e,n oc Making the fullest 

allowance for broken structural homology thus brings the observed FP closer to the 
virial plane, with the exponent of the surface brightness term consistent with the virial 
expectation. 

Key words: galaxies: fundamental parameters - galaxies: structure - galaxies: ellip¬ 
tical and lenticular, cD - galaxies: dark matter - galaxies: kinematics and dynamics. 


1 INTRODUCTION 
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It has been known for some years that elliptical galaxies are 
well represented by a two dimensional s urface in the space 
defined by their observable parameters (Brosche fc Lentei 
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1983; Djorgovski 1987a; Lynden-Bell et al. 1988a). This two 


dimensional manifold, commonly shown in the logarithmic 
3-space of radius (Re), intensity (7e), and central velocity 
dispersion (ao), has come to be known as the Fundamen¬ 


sis models, wher e metallicity ef fects result in only a small 
tilt in the plane (Worthey 1994). 


tal Plane (FP) (Djorgovski & Davis 1987; Dressier et al^ 
1987b|^ The existence of such a plane has implications for 
the formation and evolutionary processes of elliptical galax¬ 
ies. Furthermore, the small scatter about the FP makes it 
a useful tool for estimating distances to elliptical galaxies. 
Acknowledging some degree of scatter about the FP, either 
intrinsic to the galaxies and/or due to measurement errors, 
the deviati ons from the FP have been interpreted as pec uliar 
velocities (Dressier et al. 1987a; ^Dressier fc Faber 1990 ). For 
such methods of studying peculiar velocity helds and bulk 
flow motions, the exact FP used is a crucial factor. 

Under the assumptions of kinematic and structural ho¬ 
mology (i.e. identical form and shape for the kinematic and 
density profiles) amongst galaxies, the virial theorem pre¬ 
dicts that galaxies will occupy a plane in log-space amongst 
the observables radius {Re), mass surface density (rie ), and 
central velocity dispersion (tJo) (Faber et al. 1987). The 
transformation between this plane and the observed plane is 
given by the mass-to-light ratios, M/L, of the galaxies. Fur¬ 
thermore, the tight constraints on the thickness of the ob¬ 
served plane and its negligible curvature, implies that M/L 
is a power law function of the observables. The mass-to- 
l ight ratio has important implication s for galaxy formation 
( Faber et al. 1987 ; Djorgovski 1987b ) and so one would like 
to be able to test the assumptions of homology that go into 
the model. 


This has led to investigations that test the assumptions 
of kinematic and structural homology and measure the in¬ 
fluence of such a breakdown upon the FP (Ciotti, Lan zoni fc 


Renzini 1996; Hjorth & Madsen 1995). Pahre et al. (199 


also concluded that systematic departures from structural 
and dynamical homology amongst the galaxies may be re¬ 
sponsible for the departure between the observed and the 
virial form of the FP. 


Th 


Capelato, de Carvalho & Carlberg ( 1995 ), using com¬ 
puter model simulations of merged galaxies, claim that the 
slope of the FP is largely explained by broken h omolo gy in 
the velocity distri butio n. However, Ciotti et al. (1996) and 
Ciotti & Lanzoni (1996), using spherical models with the ve¬ 
locity anisotropy described by the Osipkov-Merritt formula, 
find it is unlikely that orbital anisotropy plays a major role 
in producing the FP tilt. Anisotropy in the velocity distribu¬ 
tion has also been rejected by Djorgovski & Santiago (1993) 
as the sole cause for the M/L ~ L“ relation. Ciotti et al. 
( 1996| ) suggested that a systematic trend in the shape of the 
elliptical galaxy light proHles seems to be the best explana¬ 
tion. 

This possibility has been entertained before (Djor- 
govksi, de Carvalho, & Han 1988; Djorgovski & Santiago 
1993), but differing profile shapes causing the tilt of the 
FP has not been explored with observational data. Hjorth 
& Madsen (1995) parameterised broken structural homol¬ 
ogy in their models of elliptical galaxies based on the sta¬ 
tistic al mechanics of violent relaxation (Hjorth & Madsen 


observed difference between the FP and the rela- 


1991). They found that for a given correlation between the 


tion predicted from the virial theorem is usually viewed as a 
systematic variation in the mass-to-light rati o along the FP, 
given by M/L ~ L°, where a ^ 0.25 ± 0.05 ( Faber & Jack- 


son 19'^; Faber et al. 1987; Djorgovski 1987b), and it has 


galaxy structure and its luminosity, such non-homology is 
not inconsistent with the observed FP having a constant 
M/L. Unfortunately, their data set did not enable them to 
draw firm conclusions. There is, however, a growing wave of 


1993; Pahre, Djorgovski & de Carvalho 1995). Many causes 

ogy, amongst dwarf galaxies ( 

Davies et al. 1988; 

Binggeli 

of such 

an effect have been explored and dismissed to vary- 

& Cameron 199l|), spiral bulg 

es (jCapaccioli 1987; 

Capac- 

ing ext< 

ints. Possible stellar population differences along the 

cioli 1989]; Andredakis, Peletier & Balcells 1995; Courteau, 

FP ( 

Dj 

irgovski 1987t|; pjorgovski & Santiago 1993[ [Renzini] 

de Jong & Broeils 199f:), elliptical galaxies (Michard 1985; 

& Ciotf 

i 1993) seem unlikely due to the required hne-tuning 

Caon, Capaccioli & D’Onofrio 1993; Burkert 1993; Hjorth 


of the IMF or M/L required in order to maintain the thin- 


ness of ,he observed FP. A systematic variation in stellar age 


fc Madsen 1995), and brightest cluster galaxies (Schombert 


along the plane could produce the observed tilt, but would 
require a large conspiracy in the formation of all galaxies 
of a given mass in order to be consistent with the observed 
small thickness of the FP. 

To examine stellar population effects, Pahre et al. 


198C; Graham et al. 1996). 


( 1995) tested the significance of the mass-metallicity relation 
(Guzman, Lucey & Bower 1993) along the FP, by construct¬ 


ing a near-infrared FP. Sampling the older distribution of 
stars, their K-band photometry is less sensitive than optical 
wave-bands to differences in galaxy meta llicity which appea r 
as line-blanketing effects in the optical (Faber et al. 1987). 


The near-infrared FP obtained did not differ markedly from 
the optical FP, with the deviation found being attributed 
to the reduction in metallicity effects between galaxies. This 
agrees with earlie r work by Dressier et al. (198 7b), R ecillas- 
Cruz et al. ( 199C ), and Djorgovski & Santiago (1993), where 
the metallicity effect was shown to be insufficient to explain 
the entire observed FP departure from the virial plane. This 
is also shown to be the case with current population synthe¬ 


In this paper, using a set of observed surface brightness 
profiles and measured velocity dispersions, we explore the in¬ 
fluence of broken structural and dynamical ho mology upo n 
the FP. Using the Sersic light proHle (Sersic 1968), 

the generalised form of the de Vaucouleurs law (de 

Vaucouleurs 1948, 1953) with n a free parameter, we can 
allow for structural differences between galaxies. The R^^* 
law has two parameters which act as physical scales, scaling 
the radius and the surface brightness and leaving the proHle 
shape fixed. The R^^^ profile has the additional parameter, 
n, which permits different shapes for the galaxy light prohle 
and therefore allows us to examine broken structural ho¬ 
mology amongst galaxies. The Sersic law can take the form 
of both an exponential disk (u=l) and the de Vaucouleurs 
light profile (n=4). It can also approximate a power-law and 
provides intermediate forms between these common pro- 
hle sh apes by varying its s hape p arameter n. Caon et al. 
and Graham et al. (1996), and references therein, 


(tl993| ) 


have shown that a range in galaxy prohle shapes does exist, 
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and more importantly that a systematic trend between pro¬ 
file shape and galaxy size exists such that the larger galaxies 
have le ss curvature in their light prohles than the smaller 


We obtained the effective radius, Re, and mean surface 
brightness, Ee, within Re (expressed in linear units rather 
than mag arcsec“^) by fitting an R}^^ model profile (Sersic 


galaxies. In this paper we examine the effect of this on the 1966) to the corrected luminosity profiles. The form of this 


FP. We also explore the use of volumetric galaxy parameters 
as well as using the standard projected quantities. It is pos¬ 
sible to obtain the three dimensional structural quantities 
by deprojecting the light profiles (Ciotti 1991). The 

volumetric kinematical quantity is derived from application 
of the Jeans hydrodynamical equation to the R}^^ model, 
calibrated by the projected central velocity dispersion. 

In section 2 we describe our sample of elliptical galax¬ 
ies and the derivation of their structural parameters. An 
overview of the parameters used in the construction of the 
FP and their relation to the ideal virial plane parameters is 
given in section 3. We use multivariate statistics to describe 
and construct the FP in section 4. Our results are discussed 
in section 5, and our conclusions presented in section 6 . 


profile is 
I{R) = /ecxp I —b 


(I) 


1 /n 


-1 


( 1 ) 


where le is the intensity at the radius Re- The parameter h 
is chosen so that Re is the projected radius enclosing half 
the total light from th e galaxy and is w ell approximated by 
b = 1.9992n - 0.3271 (|Capaccioli 1989|) . 

The aperture magnitude, within radius R, is obtained 
from the projected luminosity within a circular aperture of 
this radius. The luminosity within a radius R, is then cal¬ 
culated from 27r I{R)RdR. Using the change of variable 
X = b{R/Re)^^", it can be shown that 


L{R) = 2Tin—'){2n,x)IeRl 


( 2 ) 


2 DATA 


2.1 Surface Brightness Profiles 

For the purposes of this investigation we have chosen the 
s mall but very high-quality data set of Bower, Lucey & Ellis 
(1992). The aperture magnitude profiles of 26 E/SO Virgo 


cluster galaxies were selected, due to the quality of this CCD 
data (with an rms internal scatter of 0.015 mag) and as 
published dynamical data is also available for these galaxies 
(Dressier 1984, 1987; McElroy 1995). 

We took the tabulated V band profiles, and corrected 
them for the galactic absorption and redshift terms as fol- 
l ows. W e applied the K-correction as detailed by Bower et al. 
(1992), being AU=0.005 mag for Virgo. Galactic extinc tions 
in B-band (Ab) , take n from Burstein & Heiles (1984) and 
Burstein et al. ( 1987), were applied to our data such that 
Av ~ 0.75Ab (Bandage 1973). For those galaxies not listed, 
we used the NASA/IPAC Extragalactic DataBase (NED) to 
find their extinction. We also applied the corrective term for 
redshift dimming, being -0.013 mag for Virgo at a redshift 
of 0.003. So as to remove the data most affected by seeing, 
we truncated the inner profile, excluding data within 3 see¬ 
ing FWHM. Each galaxy profile consisted of 14 data points 
before allowances were made for seeing, which resulted in 
the removal of between 2 to 5 inner data points. The outer 
data p oint is given by the largest aperture fitted by Bower 
et al. ( 199^), having a diameter of 79 arcsec. 


2.2 Galaxy Parameters 


We have use d the velocity dispersion data in Table 11 of 
Bower et al. (1992), which comes from Dressier (1984, 1987). 
We note that an effective aperture of 16 x 16 arcsec was used 
to make these measurements (Dressier 1984). Four of the 
galaxies were not measured by Dressier (marked in Table 2) 
and so we have us ed th e mean central velocity dispersion 
given by McElroy (1995), correcting it t o the standard used 


by Dressier. From Table 2 of McElroy ( 1995 ), this implies 
dividing the central velocity dispersion value of McElroy by 
1.02 to obtain the value which would have been measured 
using Dressler’s system. 


where 7 is the incomplet e gam ma function. The galaxy light 
profiles of Bower et al. (1992) are in the form of aperture 
magnitudes, and we used a standard non-linear least squares 
fit (i.e. minimisation of the ordinate to define the best fitting 
line) to obtain the model parameters Re, Ate(= —2.51og7e), 
and n, from the equation 


m(r) = G — 2.5 log 7 ( 2 n, a;), (3) 

where C = —2.5 log[7ei7e27rne*’/6^"]. The residuals of the 
aperture magnitude data points about the best-fitting 77^^" 
profiles are shown in Figure There were multiple observa¬ 
tions of individual galaxies, and we show the residual profiles 
for each set of data for each galaxy. 

The model parameter Re (arcsec) was converted to kpc 
using a Virgo redshift of 0.003, and 77o = 75 km s“^ Mpc“^. 
As Re is sometimes larger than our observed outer radius 
data point of 39.5 arcsec, one cannot always simply integrate 
the profile data out to Re in order to get Ee. An alternative 
approach was developed based upon the surface brightness 
profile from the 77^^" model (see Appendix A). Such extrap¬ 
olation is valid given the assumption that the 77'^'^" model is 
suitable beyond our observed outer radial point and out to 
the model half-light radius Re- For Virgo, 39.5 arcsec corre¬ 
sponds to 2.28 kpc, and roughly 1/3 of our sample required 
an extrapolation. We averaged the model parameters when 
multiple observations existed, obtaining a single mean value 
to represent each galaxy. Agreement between different ob¬ 
servations of similar galaxies is good. This can be seen in 
Table 1, or more easily in Figure ^ where the Icry^ confi¬ 
dence intervals of the fitted model parameters overlap each 
other for multiple images of the same galaxy. 

We have plotted the logarithm of the average shape pa¬ 
rameter n, obtained from the models fitted to the galaxy 
profiles, against the logarithm of the average spatial half- 
light radius in Figure As seen in the following section, 
the spatial half-light radius is approximately 1.3477e for all 
values of n. That is not to say the spatial half-light radius is 
independent of n, as it is n that determines the profile shape 
and thus the value of Re and the spatial half-light radius. 
Also shown in the n — log Re plane are the Icr confidence re¬ 
gions for the fitted model parameters of every image. These 
are based upon the projected = 2.30 contour around 
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the best fitting solution (n, -Re, /ie), after normalising the re¬ 
duced to 1 at the fitted minimum (Figure |^. Saglia et al. 
(1996) point out that under certain conditions, such x^ con¬ 
tours can underestimate the true errors. This effect becomes 
increasingly apparent as the actual galaxy light profile de¬ 
parts further and further from the assumed model profile 
(an R^/^ law in their analysis of this effect). In fitting an 
Ri/" model, we explicitly allow for such an effect (to the 
extent that the R^^" model can correctly fit the data) and 
hence reduce the influence of this effect on our estimated er¬ 


rors. Saglia et al. ( 1996 ) note that the discrepancy is mainly 
due to large extrapolations of the fitted profile models in 
obtaining quantities like total galaxy magnitude. We do not 
compute the total galaxy magnitude, and consequently our 
model parameters are largely free from such extrapolations. 
Although Re,n is the half-light radius it does not require 
the total magnitude for its derivation, which comes directly 
from the luminosity profile fit. A few galaxies did require 
some extrapolation, but only as far as their model half-light 
radii. We note that these model half-light radii may not be 
the true galaxy half-light radii, as we don’t have detailed 
extended photometry and the R^^" model may not be ap¬ 
propriate for the faint outer parts of the galaxy, where galaxy 
envelopes or halos may result in a departure from the inner 
profile shape. The model radii obtained are however repre¬ 
sentative of the brighter inner portions of the galaxy profile 
and are useful for our analysis. Whilst there is a coupling be¬ 
tween the galaxy model parameters (although we note that 
a couple of the larger contours in Figure ^ have the effect of 
dominating this picture), such coupling cannot fully explain 
the observed trend. The value of n clearly increases with 
effective radius. _ 

Bower et al. ( 1992| ) made multiple observations of the 
26 galaxies in the data sample. Table 1 lists the model pa¬ 
rameters for every observation. The average values used in 
the construction of the FP are listed in Table 2, along with 
the observed central velocity dispersions. The values in Ta¬ 
ble 2 have been transformed from projected quantities into 
spatial quantities, as detailed in the following section. 

2.3 Possible sky-subtraction errors 

We explore the possibility that sky-subtraction errors may 
be leading to the observed departures from an p rofile. 
Even with all of the care that the Bower et al. (1992) data 
was put together, the use of small area detectors can leave 
residual sky-subtraction errors. Perhaps the larger galaxies, 
with apparently greater values of n, are a product of sys¬ 
tematically larger sky-subtraction errors which mimic de¬ 
partures from an R^^^ law. 

To investigate this, we fitted the galaxy profile data 
with the model 

_(£) -1 ( 4 ) 

which is the standard R^^^ law plus an intensity term for 
the possible sky-subtraction residual (which can be either 
positive or negative in this equation). Sky-subtraction is a 
common concern, and an incorrect value for the sky bright¬ 
ness will result in changes to the derived galaxy magnitude 
profile. The addition of an extra parameter in the profile fit¬ 
ting, namely I sky, will of course produce more accurate fits 


/(R) = 7eexp < —7.669 


as governed by the x^ statistic. The question is whether or 
not this extra parameter (actually a correction) is justifiable. 

Plotted in Figure ^ is the ratio of the x^ fit from the 
R^^^ -I- sky-correction model to the x^ fit from the simple 
R^^"^ model. Also shown is the ratio of the x^ fit from the 
r 1/" model to that from the simple model. Both of 
the additional models, containing a third parameter, clearly 
show an improvement in fit when one deviates from an n=4 
profile. The question to be addressed here is whether use of 
the sky-corrected R^^'^ model can be substantiated. We turn 
our attention to the size of the sky-subtraction correction 
that these models imply. 

The photometry com pariso n with previous authors 
work, which Bower et al. (1992) performed, puts a limit 
on the size of possible errors. Even if we attribute the en¬ 
tire disagreement between observers work to sky-subtraction 
errors, we stil l have quite stringent limits. Aaronson, Pers- 
son & Frogel (1981), using large photographic plates rather 
than small area CCD’s, obtained V-band data for 22 Virg o 
galaxies in common with the sample of Bowe r et a l. (1992). 
Rejecting 2 out lying points. Bower et al. (1992) found a 
scatter of only 0.035 mag in the comparison photometry. A 
similar scatt er of 0.036 mag was found with the comparison 
of Michard (1982) (also with 20 common galaxies). Even if 


we attribute a ll of t his scatter to sky-subtraction errors by 
Bower et al. (1992), we still have a very tight constraint 


on the allowed sky-subtraction errors in our profile models. 
Figure ^ shows that the majority of sky-subtraction errors 
are notably larger than permitted by our quite liberal error 
allowance. So while the addition of a sky-subtraction correc¬ 
tion produces improved fits to the data, it is not justified. 

The introduction of this additional parameter has the 
success it does because it is mimicking, to some degree, the 
effect of some other real physical process or situation. Whilst 
an R^^'^ law is being fitted in these models, it is being done 
coincidentally with a term that effectively re-shapes the pro¬ 
file to obtain an R^^* profile. In this fashion, the model in¬ 
corporating a sky-correction term is able to allow for real 
differences in galaxy profile shapes, as better represented 
by the R^^" law (Caon et al. 1993; Graham et al. 1996) 
In Figure we show that the R^^^ + sky-correction model 
is only superior to the R^^" model for unacceptably large 
under-estimations of the sky background level. Also, 13 of 
the points in Figure ^ for which the R^^'^ + sky-correction 
model is superior to the R^^" model belong to just one 
galaxy (NGC 4552), which was observed 14 times; 6 of these 
points occupy the cluster of 8 points that favour the R^^’* 
-I- sky-correction model. Removal of NGG 4552 eliminates 
the impression that use of an R^^'^ + sky-correction model 
is generally better than the R^^" model when one is dealing 
with the larger galaxies (n>4). 

In further support of our conclusion that sky- 
subtraction errors are not the effect responsible for the range 
in profile shapes seen, Capelato et al. (1995) also found a 
correlation between n and Re. The R^^’’ models they fitted 
to their computer simulated elliptical galaxy merger rem¬ 
nants, have no uncertainty as to sky levels, thus eliminating 
this possibility. There is however another possibility as to 
the cause of the range in our profile shapes. SO galaxies are 
known to possess a disk as-well as their central bulge. If our 
sample of galaxies contained significant disks and a range 
of bulge-to-disk ratios, then this might explain the apparent 
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broken structural homology (Saglia et al. 1996). Whilst this 
could be so, an argument against such an origin is that the 
galaxy models of Capelato et al. (1995) were disk free, as 
was t he wo rk with Brightest Cluster Galaxies by Graham 
et al. (1996), where a large range in profile shapes was still 
evident. Also, given the observed trend with n and Re, one 
would have to claim that the disk-to-bu lge ratio increased for 
the big ger galax ies, which it does not (Andredakis, Peletiei 
& Balcells 1995). 


3 THEORY 


3.1 Overview 

A good summary of the relation between the observed and 
the v irial p lanes can be found in Djorgovksi, de Carvalho & 
Han (1988), and we adopt their method of analysis in what 


follows. From the virial theorem, the connection between the 
potential and kinetic energy can be written as 


GM 
< R> 


= k, 


<V^ > 


(5) 


with M and <R> being a measure of the galaxy mass and 
radius, respectively, such that GM/<R> is the galaxy po¬ 
tential energy and <F^>/2 is the mean kinetic energy per 
unit mass for the galaxy, ke is a virialisation constant re¬ 
flecting the degree to which the system is virialised. Galaxy 
mergers and the process of violent relaxation redistribute the 
orbital energy of the stars. Galaxies may exist in a quasi¬ 
static equilibrium, with varying degrees of energy dissipation 
or in varying stages of relaxation, and hence varying values 
of ke may exist before a galaxy settles d own into virial equi¬ 
librium (fce=2) (Djorgovski et al. 1988). In addition to this, 
the ke term also allows for possible neglect of rotational en¬ 
ergy: in our use of central velocity dispersions to measure the 
kinetic term, <V^>, we have ignored possible contributions 
from the rotational energy of each galaxy. 

A choice now arises regarding which observables one 
uses to fit for the quantities M, <R> and <V^>. 


3.2 Structure 

For elliptical galaxies, an observed radius R (such as Re) is 
typically used to approximate <R> such that 

R = ku < R > . (6) 

The parameter k^ allows for the possibility that our ob¬ 
served projected radius R may not be equal to the physi¬ 
cal radius we wish to use. kpt depends on the mass-density 
structure of the galaxy, allowing for variations in the distri¬ 
bution of matter from one galaxy to the next. It is possible 
that different galaxies formed with a different mass-density 
structure. It is also likely that evolutionary processes such 
as galaxy merging and tidal interactions have further modi¬ 
fied the fcfl term. By addressing the breakdown of structural 
homology between galaxies, we can explore the variations in 
the kfi term. These variations, if ignored, introduce system¬ 
atic errors into the parameterisation of elliptical galaxies, 
and hence into the construction of the FP. 

The central velocity dispersion, no, is used to represent 
the kinetic energy term such that 


no" = kv<V^ >, (7) 

with kv reflecting the dynamical structure within the 
galaxy. Like the ka term, the kv term reflects the forma¬ 
tive and various evolutionary processes that have shaped 
the galaxy. It incorporates the anisotropies of the velocity 
dispersion tensor, including the differing degrees of random 
and circular motion in the galaxy, and differing velocity dis¬ 
persion prohles. 

Lastly, we have 

(M/L)T.R^ = (8) 

kL 

where (M/L) is the global galaxy mass-to-light ratio, and 
kL incorporates the luminosity structure within the galaxy. 
E is an intensity term, obtained from the htted light profile. 
Substitution into equation ^ gives 

R^KsroIt.-\M/L)-\ (9) 


with 


Ksr = 


he 


2Gkfiki,kv 


( 10 ) 


being the combined structural term appr opria te for this 
expression, as given by Djorgovski et al. (1988). Without 


knowledge of the Ksr term, it has, in the past, been assumed 
a constant, and the observed departures of the FP from the 
relation R oc have been attributed to (M/L), per¬ 

haps with an additional comment as to uncertainties about 
assumed structural and dynamical ho molo gy. This paper ad¬ 
vances the work of Djorgovski et al. (1988), where they pro¬ 


ceeded under the assumption of homology (i.e. that all the 
k’s are constant) - we do not. Through fitting an pro¬ 
file which specifically allows for structural differences, there 
is no assumption of structural homology. We point out that 
whilst the R^^" profile covers an extensive range of proHle 
shapes, it does not cover every possible profile shape as that 
would require an endless number of parameters in the ht¬ 
ted prohle model. Thus, we only explore broken structural 
homology as far as the model allows, which turns out 

to be quite sufhcient. Allowing for variability in the prohle 
shapes, we can now better address the ku and fci, terms, and 
also have a better representation for the ky term. 

We explore the use of spatial (i.e. three-dimensional) 
quantities to represent <R>,M, and <V^>. We take the 
spatial half-mass radius, ri/ 2 , which is equivalent to the spa¬ 
tial half-light radius for a constant M/L, to represent <R>. 
We compute the total galaxy luminosity, and use the as¬ 
sumption of constant M/L to obtain half the total galaxy 
mass which is used to represent M. The spatial velocity dis¬ 
persion at ri /2 is acquired using the Jeans equation (see 
Section 3.2) and is used as our estimate of <V^>. These 
quantities are not directly observable, but they can be com¬ 
puted from other observable quantities given a knowledge of 
the k-terms. 


Giotti (1991) related the luminosity form of the Sersic 


law to its underlying mass-densi ty structural form, by solv - 
ing the Abel integral equation (Binney & Tremaine 1987) 


which relates the projected/deprojected expressions. He also 
computed the spatial half-mass radius ri /2 for differing ex¬ 
ponents n, and found that this radius, divided by the pro¬ 
jected effective radius Re, is practically a constant value 
(ranging from 1.339 to 1.355 for n=2 to 10). We note a sim- 
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ilar figure was also f ound for the Jaffe (1983) models (1.311) 
and the Hernquist (199C) models (1.33). Thus, the spatial 
half-light radius (in units of Re) is largely independent of 
the variant forms of the mass-density distributions used to 
describe the galaxy and so ka can be taken to be a constant 
with 


R = iZe,n « (1/1-34) <R>= 0 . 75 ri/ 2 . (11) 

This is not to say that ri /2 is independent of n and there¬ 
fore does not depend on the galaxy profile shape. What it 
does say is that ri /2 is dependent on the projected half-light 
radius, which is very much dependent on the profile shape 
and hence the value of n. 

The structural term /cl is not as simple, but it can be 
computed for different values of n. Using M = Mtot/2 and 
{M/L) = {Mtot/Ltot) in equation |, upon substitution of Ltot 
from equation ^ (where the incomplete gamma function, 7 , 
becomes the gamma function T), we have 

kL ~ 7rnr(2n)e*’/6^". (12) 

As b is well approximated by a simple linear function of n 
(l<n< 10 ), we have an expression for the luminosity struc¬ 
tural term in terms of the shape parameter n. Allowances 
for different galaxy luminosity structures are made to the 
mean intensity term E, which is fed into the FP as the 
surface brightness parameter. This is achieved by using 
fie — 2.51og(A:L), rather than /ie (the surface brightness at 
Re) which does not fully represent the luminosity structure 
(as detailed in equation ^). If all galaxies did obey an R^^^ 
law, then the difference between these terms would just be a 
constant and use of one term over another would not effect 
the tilt of the FP. 

The mean surface brightness, (/r)e, within the effective 
radius, Re, is related to the intensity by —2.51og(/)e, where 
a)e = Ltot/2'KRe from Appendix A. In Appendix A, we also 
see that (/r)e — [^e — 2.51og(fci,)] = 2.51og(7r), independent 
of the value of n. Thus (/r)e can be used in the construction 
of the FP, instead of fie — 2.51og(fci:,). The intercept of the 
FP will be effected but not the slope of the plane. 

Thus, to the extent that the R^^'^ model represents the 
observed variations in structural form, we are able to ad¬ 
dress two of the components that go into making the Ksr 
term which determines differences between the observed and 
the virial plane. The model parameters. Re and {fi)e 

implicitly allow for these variations in profile shape which is 
explicitly represented by the third model parameter n. 


3.3 Dynamics 


As we do not confront the issue of orbital anisotropy, we 
cannot fully address the kv term, although we do partially 
deal with it by tackling the problem of differing velocity dis¬ 
persion profile shapes due to differing luminosity profiles, 
as represented by the R}^^ model. Dynamical homology im¬ 
plies that each galaxy has the same fall off in velocity disper¬ 
sion away from its center. It is therefore obvious that broken 
structural homology necessitates broken dynamical homol¬ 
ogy. Working with the assumptions o f sphe rical symmetry 
and isotropic pressure support, Ciotti ( 1991 ) has shown how 
both the spatial velocity dispersion profile and the projected 
line of sight velocity dispersion profile change with radius 


for a series of R^^'^ models with a range of shape param¬ 
eters from 2 to 10. We do not deal with possible dynami¬ 
cal anisotropies within each galaxy, as this requires detailed 
spectral mapping and/or assumptions about the ratio of the 
radial to azimuthal velocity dispersion (Ciotti, Lanzoni & 
Renzini 1996). 

There is of course the additional concern of an aper¬ 
ture effect. Even if all galaxies did have the same velocity 
dispersion profiles, as one samples galaxies that are further 
away, the fixed aperture size on the sky will sample a larger 
portion of the galaxy. T his ha s been addressed by Dressier 
(1984) and Davies et al. (1987), and is implicitly taken into 
account in what follows. However, the problem is more com¬ 
plicated than this. What one would actually like to do, is not 
be consistent in using the same physical aperture size on 
each galaxy, but to sample each galaxy with the same ratio 
of aperture size to half-light radius^ (Wegner et al. 1996). 
That is, to always measure the velocity dispersion within 
say, IRe or 0.2i?e or some fixed ratio of Re- It is clear that 
use of core velocity dispersion measurements is inadequate 
in representing a galaxies overall ki nemat ic state. 

Extending the work of Ciotti (1991) we compute the 
projected aperture velocity dispersion as a function of the 
ratio of the aperture size to the projected half-light radius, 
as shown in Figure (The details are left for Appendix B). 
From these models, the ratio of the projected aperture ve¬ 
locity dispersion (of aperture size equal to that used to ob¬ 
tain the central velocity dispersion at the telescope) to the 
spatial velocity dispersion at ri/ 2 , can be used to correct 
the measured central velocity dispersion to the spatial half¬ 
mass value. This ratio is shown in Figure ^ Depending on 
the aperture size used (actually the ratio of the aperture size 
to the projected half-light radius), use of the measured cen¬ 
tral velocity dispersion by itself may result in departures of a 
factor of 2 or 3 from the spatial half-mass velocity dispersion 
measure. 

It is possible to estimate the galaxy masses, using our 
measures of galaxy size and velocity dispersion. Once we 
have fitted the galaxy luminosiW profile with an law, 

the only unknown in equation B 8 is {MjL), where we have 
used our measured central velocity dispersion for Uap,n^ap) 
and the value predicted by theory, shown in Figure for 
'^ap,n(^ap). Substituting in the total galaxy luminosity for 
L, from equation ^we can then solve for the total galaxy 
mass Mtot- Figure y shows the total galaxy mass obtained 
in this way, plotted against the shape parameter n. 

We also use the model velocity dispersion profiles and 
the observed central velocity dispersion to compute the aper¬ 
ture velocity dispersion extended to infinity, atot,n- Ciotti 
(1994) shows this quantity to be equal to one third of the 
virial velocity dispersion (in the case of spherical systems) 
and independent of any possible orbital anisotropy. 

Finally, we remind the reader that thi^aper does not 
treat the fourth parameter ke, in equation M which allows 
for different degrees of virialisation and the neglect of rota¬ 
tion. 


t This correct approach would imply a differe nt conversion of the 
central v elncitv dispersion in McElroy ( |l995 ) to the standard of 
Dressier (1984). 
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4 FUNDAMENTAL PLANE CONSTRUCTION 


With the FP relation described as 


We construct four planes. The first one is the standard plane, 
obtained by fitting profiles, and uses the projected ef¬ 
fective radius i?e, the mean surface brightness within this 
radius, Ee, and the measured central aperture velocity dis¬ 
persion (To- This plane is used as our reference plane, to 
see how the planes produced by alternate choices of vari¬ 
ables change. The second plane that we construct allows for 
differences in galaxy structure. We fit an R}^'^ profile, and 
use the spatial half-light radius, the associated mean surface 
brightness within this radius, and the central aperture veloc¬ 
ity dispersion. Such a plane will tell us the extent to which 
use of an R^^'^ law for all galaxies is responsible for the tilt of 
the FP. The third plane treats broken structural homology 
and also allows for the fact that a fixed aperture measure 
for the velocity dispersion will sample different ratios of the 
area enclosed by Re, as Re itself varies between galaxies. 
Here we transform the central aperture velocity dispersion 
measure into the spatial half-light velocity dispersion cri/ 2 ,n- 
The fourth plane we construct uses the spatial structural 
parameters from the R^^" model and an infinite aperture 
velocity dispersion term. 

A Principal Component Analysis (PCA) was performed 
on each 3 dimensional data set, (log no, log Re, 4 , Se, 4 ), 
(logCro,logRe,n,Ee,n), (log , log Re,n, Ee,n ) , and 

( log (J tot.n, log Re.n, Ee.n). The code from Murtagh & Heck 
(1987) was used to show the degree to which the data are de¬ 
fined by a plane in our 3-space of observables. Table 3 shows 
that greater than 98 per cent of the variance in the data lies 
in a plane, confirming the appropriateness of constructing a 
plane to describe the properties of elliptical galaxies. 

The method of construction of the plane is not an obvi¬ 
ous task and many different approaches have been tried in 


log Re{kpc) = A log cr — 0.4REe -I- C, 


(13) 


the past (Djorgovski fc Davis 1987; Lynden-Bell et al. 1988b 


Lucey, Bower fc Ellis 1991a| ). In the fitting of a plane to 3- 
space data, one has considerable choice in exactly what it is 
that one minimizes in order to obtain the ‘best’ plane. One 
gets a good idea of the p roblem from Isobe et al. (1990) and 


Feigelson & Babu ( 1992 ), where 6 possible approaches to the 
fitting of a straight line in a 2-space data set is detailed. The 
situation is of course worse for higher dimensions. However, 
extending the ideas learned in 2-space to 3-space, we have a 
better idea of how to approach the problem. 

Our current work wishes to establish the FP for the 
sake of comparison with what theory, i.e. the virial the¬ 
orem, predicts. Given this, we do not distinguish any de¬ 
pendent/independent variables, and so we cannot build our 
plane by minimising the residuals in one variable alone. An 
approach that considers the residuals in all variables shall 
be needed, and from the 2 dimensional case studies of Isobe 


et al. (1990) and Feigelson & Babu (1992), we see that the 
bisector method of regression is most suited for our needs. 
For studies requiring symmetrical treatment of the variables, 
the bisector method of linear regression was found to be the 
preferred method, having a greater certainty on the slope of 
the fitted line than other methods, such as orthogonal re¬ 
gression. In this technique, we construct the plane 3 times, 
each time minimising the residuals of a different observable 
against the other two observables. The mean slope angles, 
obtained from the average of the three planes, is then used 
to represent our ‘best’ FP. 


the virial theorem predicts that A should be 2 and B should 
be -1. The observed FP, as determined by the bisector 
method, for our four data sets described above, is shown in 
Table 4, with the associated uncertainty being the standard 
error. 

Because there is no natural ‘best’ method of construc¬ 
tion for the FP, we used an additional method, pioneered 
by Djorgovski & Davis (1987), as a check. First, the ‘mix¬ 
ing’ value b that gives the highest linear correlation coeffi¬ 
cient between log Re and (log cr-|-feEe) is obtained. We used 
the bootstrap method to estimate the standard error on b. 
It is then possible to treat the problem as a 2-dimensional 
least-squares proble m. The program SLOPES, from Feigel¬ 
son & Babu ( 1992| ), was used to compute the best fitting 
line for a multitude of methods: ordinary least squares of 
y on X OLS(y|x), and the inverse OLS(x|y), the line which 
bisects these two, minimisation of the perpendicular residu¬ 
als (known as orthogonal regression) the geometric mean of 
OLS(y|x) and OLS(x|y) (referred to as the reduced major 
axis regression line) and finally the arithmetic mean of the 
OLS lines. The reader is referred to their paper for a discus¬ 
sion of each method, for which the slopes and errors are com¬ 
puted using three techniques. The first technique uses the 
linear regression formulae, but as these can underestimate 
the slope variance for small data samples, a bootstrapping 
simulation and a jackknife analysis are also performed. The 
various solutions are displayed in Table 5. The value of A in 
equation is explicitly shown and the value B is calculated 
such that B = —2.5 x b x A, where b is the ‘mixing’ value 
mentioned above. 

To represent the method of analysis using the mixing 
parameter, we use the orthogonal regression solutio n from 
Table 5, consistent with the method of Pahre et al. (1995), 
as the one method from the possible 6 . The bisector method 
should perh aps be used, given the discussion in Feigelson & 
Babu (1992), but the differences between these two planes 
are within the errors. This orthogonal method gives sim¬ 
ilar planes to those found with our 3D method of analy¬ 
sis. Figure shows the various FPs constructed using our 
3D bisector method. Figure uses the standard data set 
from the R^^'^ fits and the central aperture velocity disper¬ 
sion measures, whilst Figure [lo| b uses the data set from the 
R^^" fits and the central aperture velocity dispersion. Fig¬ 
ures pT| c, pT| d use the spatial parameters from the R^^" fitted 
model along with the spatial half-light velocity dispersion 
term and the infinite aperture velocity dispersion term re¬ 
spectively. 


5 DISCUSSION 

The additional parameter, n, in the R^^” light profile mod¬ 
els, is shown to have a physical association with the galaxies’ 
other properties, such as size and mass, as seen in Figure ^ 
and Figure The bigger galaxies have larger values of n, 
and thus less curvatur e in t heir light profiles, as prev iously 
shown by Caon et al. (1993) and Graham et al. (1996). Not 
allowing for variations in profile shape can have a substantial 
systematic effect on the derived half-light radius and asso- 
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dated surface brightness term. Figure [ll] shows the ratio of 
the half-light radius derived using the model to that 
obtained using the model, plotted against the model 

value of n. One clearly sees that assuming an law for the 
larger galaxies will result in an underestimate of the effec¬ 
tive radii. So too, use of the law on the smaller galaxies 
will lead to over-estimates for the galaxies’ effective radii. If 
the galaxy luminosity profile is best described by a value of 
n greater than 4, the effective surface brightness will be too 
bright if an profile is applied. Likewise, the effective 
surface brightness term from an profile is fainter than 

that from the R^^'^ profile when n is less than 4. Allowing for 
structural differences has the result of producing a greater 
range in both radius and surface brightness. 

With the introduction of this additional parameter, we 
performed another Principal Component Analysis to de¬ 
termine the optimal number of dimensions in the 4-space 
(n, log i?e,n, Se,n, log (To). The major eigenvector contained 
67 per cent of the variance, and the second dominant eigen¬ 
vector had 27 per cent. Thus, our 4-space of observables is 
still well represented by a two dimensional plane, encom¬ 
passing 94 per cent of the total variance (use of log(Ji/ 2 ,„ 
and logcTtot.n gave similar results). Therefore, a hyperplane 
does not seem warranted. This is not surprising due to the 
strong correlation seen between n and log Re in Figure p. 

If one was to use the FP as a distance indicator, then 
its method of construction would be different to the method 
used here. In such a case, one would wish to predict the 
value of one distance-dependent variable from the measure¬ 
ment of other distance-independent variable(s), and so one 
should use OLS(y|x), where y is the variable to be predicted. 
With this OLS(y|x) fit we obtained an rms scatter about the 
plane of 0.09 in log Re when an 7?^^^ law was applied with 
(Jo, and an rms scatter of 0.10 when an 7?^^" law and spatial 
structural quantities were used with (Jtot,n. The uncertainty 
on these rms scatters is around 0.02, so there is no signif¬ 
icant difference in the scatter about the FP using the two 
methods. We also measured the intrinsic thickness of the FP 
by computing the rms orthogonal offset of the data points 
about the plane; here we assume no dependent/independent 
variables. These rms values are shown in Table 4. 


The virial theorem predicts the relation R (x (j^E“^, 
when the surface brightness, E, is expressed in linear units 
rather than mag arcsec“^. Observed planes for cluster 
galaxy data, using 7? = 7Ze,E = Ee, and (t = (tq (the pro¬ 
jected central velocity dispersion) have been in good agree¬ 
ment with each other in the past. Dressier et al. (1987b) 


found R oc a 


1.33v^-0.83 


for 6 combined galaxy clusters. 


with Lucey et al. (1991a) obtaining R oc (j^'^^E“'^’^® for a 
Virgo s ubsample (/re(B)<22.00) of 16 ellipticals. Lucey et al. 
( 1991b| ) found a similar result for their V-band study of 51 
ellipticals in Coma, with R oc (ji-23 j^-o .82 ^ larger study 

of 11 clusters by Jprgensen, Franx & Kjaergaard (1996), 
taken in Gunn r, produced an identical result. 


A recent near-infrared study of 5 near by clusters 
(59 elliptical galaxies) by Pahre et al. (1995) gave R oc 
^i.44j^-o.79 steeper slope of the plane was attributed to 

a reduction of metallicity effects in the near-infrared band¬ 
pass compared to the optical. Intriguingly, they obtained 
the notably different R oc for their Virgo sam¬ 

ple of galaxies. As this sample only numbered 8 galaxies. 


this result does not carry much weight on its own, but does 
warrant further investigation. 

As shown in Table 4, allowing for structural homol¬ 
ogy through use of an 7?^/" profile, and maintaining use 
of the central velocity dispersion, actually increases the de¬ 
parture of the FP from the virial plane. Whilst the radial 
term changes in the sense that it reduces the departure, the 
surface brightness term works to increase the departure and 
proves to be the dominant effect. However, changes to the 
structure of a galaxy imply changes to the velocity field of 
the galaxy. These changes can be obtained through applica¬ 
tion of the Jeans equation for a pressure-supported system. 

One problem with the use of central velocity disper¬ 
sions is that they sample different portions of the half-light 
radius for each galaxy. To be consistent in our measuring of 
galaxy parameters, we computed all parameters at the spa¬ 
tial half-light radius. This meant converting central aper¬ 
ture velocity dispersions into the associated value at the 
spatial half-light radius (see Appendix B for details). Do¬ 
ing this, and using the spatial parameters from the 7?^^" 
profiles, we found that the FP constructed from these quan¬ 
tities more closely resembled the virial plane. The exponent 
on the surface brightness term in this expression for the FP 
is consistent with the virial theorem prediction of -1. Past 
differences between this exponent and the virial expectation 
can be explained by systematic departures from structural 
and dynamical homology and the use of projected rather 
than volumetric velocity dispersion. The exponent on the cr 
term is closer to the value expected from the virial theorem, 
changing from 1.33 to 1.48 as we allow for such changes, but 
it does not reach the expected value 2. 

Although use of the spatial half-light velocity dispersion 
for all galaxies is consistent and methodical, and preferred to 
the use of central velocity dispersion measures, it is still arbi¬ 
trary. We have therefore derived the total aperture velocity 
dispersion (i.e. the value of the aperture velocity dispersion 
at infinite radius), which for spherical systems is equal to 
one-third of the (three dimensional) virial velocity dispersion 
(Ciotti 1994). The FP that allows for differences in galaxy 
profile shape and uses the total aperture velocity dispersion 
more closely resembles the virial plane, R oc (j^E“^, than 
the standard observed plane obtained with R^^^ law param¬ 
eters and central aperture velocity dispersions. The respec¬ 
tive planes are described by 77 (x fj^-^‘^tz0.ii^-o.93±o.o8 

Ro,^1.33±0.10^-0.79±0.11 

Capelato et al. (1995) found, using computer simula¬ 
tions to build and merge elliptical galaxies, that as they 
increased the aperture size in which they measured the ve¬ 
locity dispersion term, the difference between the FP they 
built and that predicted by the virial theorem was reduced. 
Although the initial conditions of their models were some¬ 
what limited and their range of merged galaxies does not 
fully represent the range of elliptical galaxies we see today, 
their work strongly indicates that the assumption of homol¬ 
ogy amongst elliptical galaxies is invalid. They effectively 
held constant the stellar population, did not alter the IMF 
or M/L ratios, and have no metallicity effects to worry about 
in their dissipationless models. The tilt of the FP in their 
models is therefore not due to any of these effects. Their 
work implied that most of the nonhomology lies close to 
the center of their merged galaxy models, and hence use of 
larger apertures gives a more accurate picture of the viri- 
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alised system — thus the reduction between the observed 
and virialised planes found with the use of a larger aper¬ 
ture. A similar physical effect could be occurring with our 
use of an effectively larger aperture, albeit an infinite aper¬ 
ture dispersion based on transforming the central aperture 
dispersion using the light models, compared to the av¬ 
erage aperture size of only 0.2Re that has been used in the 
past. 

Although our model velocity dispersion profiles allow 
for different distributions, as described by the model, 
they are built upon certain assumptions, and we have not 
dealt with possible orbital anisotropies, as would have been 
done implicitly in the computational work of Capelato et 
al. (1995). In examining the structural and dynamical devi¬ 
ations from homology between the progenitor and merged 
galaxy models, Capelato et al. (1995) claim the nonhomolo- 
gous velocity dispersion profiles as the main cause of the tilt 
to the FP. However, their merger remnants only spanned a 
range in n from 1.2 to 3.4, so they have not really explored 
the full range of structures observed in nature, and cannot 
entirel y exclude this as a contributing factor to the tilt of 


the FP. 

Ciotti et al. (1996) and Ciotti & Lanzoni ( 1996 ) con¬ 
cluded that velocity anisotropy cannot cause the observed 
FP tilt by itself (i.e. if structural homology is assumed). 
This result was derived from a range of theoretical galaxy 
models (including the model) and introducing vary¬ 

ing amounts of orbital anisotrop y as described by the 
Osi pkov -Merritt parameter isatio n (Osipkov 1979; Merritt 
1985 |). I Ciotti & Lanzoni (1996) did however suggest an 
observational test for possible velocity anisotropy effects, 
namely the use of large apertures to measure the velocity 
dispersion term, an idea supported by Capelato et al. (1995) 
and inc orporated in this study. Our work suggests that the 


combin rtion of broken homology in both the surface bright- 
ness and dispersion profiles of elliptical systems is partly 
responsible for the departure of the observed FP from the 
predicted plane. 

These results have implications for most studies that 
have been done to date using the FP, or the D — a relation, 
because such studies have assumed structural homology and 
used projected quantities. Studies of the peculiar motions of 
ellipti cal galaxies in the local Universe_j|L^ndemBelle^_M 
1988b 


Lucey & Carter 


Dressier & Faber 1990) now 


warrant re-examination, as such measurements may be sys¬ 
tematically affected by biases in the distance indicator rela¬ 
tions. Furthermore, studies of environmental dependence on 
the D — a rel ation and d ifferences between field and clus- 
ter elliptic als ( SiF^MT89| ^jorgovski^^e1^^aT]T88| Lucgy^ 
al. 199] ^ Burstein, Faber & Dressier 1990; de Carvalho & 
Djorgo'^ ski 1992) can be re-examined. If the galaxy struc¬ 
ture is dependent on the galaxy environment, then the D — a 
relation may also be dependent on environment. However, 
in fitting an profile, we have already taken into ac¬ 

count different galaxy structures (and so possibly different 
galaxy environments). Thus the new FP which we derive 
may implicitly take into consideration some environmental 
effects, although not stellar population differences (Guzman 
& Lucey 1993; Guzman et al. 1993) or intrinsic ally differ¬ 
ent M/L ratios in different clusters ( Kaiser 198^ . It will of 
course be of interest to obtain data from other clusters, and 
field galaxies, to investigate such possibilities. 


Working with the simplification that the Ksr term in 
equation ^ is now a constant (although it still contains the 
fee and kv terms which we have not properly dealt with) we 
can compute the mass-to-light ratio, assuming it is a power 
law function of the observables. For the plane constructed 
from the law model and using the central aperture ve¬ 
locity dispersion measure, we have that M/L ^ ^ 0.25 j^-o. 06 . 
whereas for the plane that was constructed allowing for bro¬ 
ken structural homology and using the volumetric quanti¬ 


ties, we find that M/L ~ 


or that M/L ~(j 


1/2 


For the plane that used model parameters and the 

total aperture velocity dispersion measure we have that 


M/L 


\ Of course it may be that M/L is con¬ 


stant and it is the Ksr term in equation ^ that varies with 
L and E or cr, as mentioned previously, or it may be a that 
terms are varying. 

The tilt of the FP is probably due to a com bination 
of factors. For instance, whilst Pahre et al. ( 1995| ) did not 
recover the virial plane in their construction of the near- 
infrared FP, th ey did measure a change in a, with M/L ^ 
M°', from 0.23 (Lucey et al. 1991b; de Carvalho & Djorgov- 


ski 1992) in the optical (V-band) to a=0.16 (K-band). They 


found that the exponent on the velocity dispersion term. A, 
increased by 0.19T0.06 dex from the V-ba nd to the K-band 
(cf. -l-0.29T0.il, Djorgovski & Santiago (1993)). 


Whilst we have dealt with the effects of seeing on the 
galaxy light profile, by excluding data within the central ra¬ 
dius of 3 FWHM, we have not dealt with the seeing effect on 
the central velocity dispersion measurements. As an effective 
16x16 arcsec aperture was used to obtain the velocity dis¬ 
persion for our sample (Dressier 1984, 1987), seeing should 
not be a problem. However, studies using a single smaller 
aperture (say 2x4 arcsec) may be effected by seei ng if there 
is a steep velocity dispersion gradient present (Whitmore 


198C; Jarvis & Peletier 1991). 


Several other issues must be kept in mind: the galax¬ 
ies may not be virialised; we have not dealt with the mat¬ 
ter of rotational support (as opposed to pressure support); 
the issue of broken dynamical homology has not been fully 
treated, as we assumed orbital isotropy; stellar population 
variations with galactic radius may result in (M/L) varying 
with radius; or light may not trace mass in a simple fashion. 
Any combination of these effects may work in symphony to 
account for the FP tilt. 


6 CONCLUSIONS 

A range in profile shapes is found to exist for the sample of 
26 elliptical galaxies used in this study. The shape param¬ 
eter n is found to have a physical association with galaxy 
size and mass, in the sense that a larger value of n corre¬ 
sponds to a larger galaxy. Allowances for different profile 
shapes and the use of infinite aperture velocity dispersion 
measurements, which allow for different velocity dispersion 
profile shapes, results in a modification to the FP. While 
more clusters need be studied, our work shows that broken 
structural homology is not solely responsible for the tilt of 
the FP, but allowing for variations caused by such structural 
differences does reduce the tilt. It is likely that the departure 
of the observed FP from the virial plane will be further re¬ 
duced by using infrared luminosity profile measurements to 
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minimise stellar population effects. Our best estimate of the 
FP using the parameters from an profile and the total 
aperture velocity dispersion is R (x 

The exponent of the surface brightness term is consistent 
with the prediction of the virial theorem. 
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Figure 1. The residual V-band aperture magnitude data (beyond 
3 FWHM, and taken from Bower, Lucey & Ellis (1992)) about the 
best fitting and profile models is shown, being open 

triangles and solid squares respectively. 

Figure 2. The average shape parameter n, from the models, 

is plotted against the average spatial half-mass radius for the 26 
E/SO Virgo galaxies from Bower, Lucey, Ellis (1992). 

Figure 3. The shape parameter n, from the models, is 

plotted against the projected effective radius for all the the images 
(including multiples). Also shown are the Icr error ellipses, as 
described in the text. 

Figure 4. The value of the best fitting R^^'^ (filled square), 
and {R^^^-\- sky-correction) (star), model, normalised against the 
X^ value of the best fitting model, is plotted against n. The 

ratio is of course close to 1 near n = 4, and is smaller when away 
from n=4, indicating a superior fit of the model to the data. 

Figure 5. The x^ value of the best fitting {R^^^-\- sky-correction) 
model, normalised against the x^ value of the best fitting R^^^ 
model, is plotted against the size of the sky correction expressed 
as the associated change in the half-light magnitude. 

Figure 6. The ratio of the x^ fit from the best fitting R^^'^ model 
to the best fitting i?t/4 model with a sky-correction, is plotted 
against the size of the sky correction expressed as the associated 
change in the half-light magnitude. 

Figure 7. The spatial (cr(r), solid), line-of-sight {cris{R), dashed), 
and aperture (crap(i?), dotted) velocity dispersion profiles are 
shown for a range of shape parameters n for the luminosity 

profile. The abscissa is such that s = r/Re and r} = R/Re. 

Figure 8. The ratio of the aperture velocity dispersion, as calcu¬ 
lated in Appendix B, to the spatial half-mass velocity dispersion, 
as a function of projected aperture size. The abscissa is such that 
7] = R/Re 

Figure 9. Total galaxy mass plotted against the shape parameter 

n. 


Figure 10. a) The FP constructed using an law to de¬ 

rive the spatial half-mass radius (R) and mean surface brightness 
(/i); and using the central projected aperture velocity dispersion 
measure for a. b) Similar to (a) except that profile parame¬ 

ters were used, c) Similar to (b) except that the spatial half-light 
velocity dispersion term is used instead of the central aperture 
velocity dispersion, d) Constructed using the parameters from 
an law and the total (infinite) aperture velocity dispersion 

term. 

Figure 11. Ratio of the derived effective radii from the R^^'^ and 
the model is shown against the shape parameter n. A clear 

systematic effect of under and over estimation of the effective 
radius is evident when applying the R^^^ law to large and small 
galaxies respectively. 


© 0000 RAS, MNRAS 000, 000-000 



Broken Homology and the Fundamental Plane 


13 


Figure Al. The structural luminosity term Kl (s.t. - 

2.51og(iCx^) = ({i)e — fie — 2.5 log(7r)) is plotted as a function of n, 
for Sersic’s law. The dashed line aids the eye for the case 

where n = 4, giving the result for the de Vaucouleurs law. 


APPENDIX A: DERIVATION OF (^)e FROM 
AN CURVE OF GROWTH 

From equation 1^ the Sersic law (3ersic 1968), describing the 
projected light distribution for elliptical galaxies, bulges and 
disks is given by 




(Al) 


The average intensity, <7>e, within the effective radius, Re, 
of the model is obtained from 

f IdA he’’ 

This expression can be simplified by letting x = h{R/Re)^^", 
giving 

<I>e=Ief{n), (A3) 

where 


f{n) = 


2ne‘’ 

b2n 


— X 2n — l 7 

e X ax. 


(A4) 


Now as b was chosen such that Re is the radius containing 
half of the total light from the galaxy, we have 


fin) = 


62 " 


—X 2n—l 

e X 


= V’ 


whe re fc z, is the luminosity structure term mentioned in sec¬ 
tion 3.2. By taking the logarithm, we obtain 


(p)e = /ie - 2.51og[/(n)]. (A6) 

The value of fie comes straight from the fitted R^^" model, 
and f{n) is computable and shown in Figure Al for a range 
of n that covers typical observational data. This r esult is in 
agreement with the previous work of Caon et al. (1994). 
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APPENDIX B: VELOCITY DISPERSION 
PROFILES 

The spatial luminosity density v(r) can be computed from 
the projected luminosity profile I{R) of the galaxy, such that 


i/(r) = 


-1 


dI{R) dR 
dR 


(Bl) 


where R is th e projected radius and r i s the spatial radius 
of the galaxy (Binney fc Tremaine 1987). 


This has been done for the R^ ^'^ profile (Poveda, Itur- 
riaga & Orozoco 1960; Young 1976), and the generalisation 
to the profile (Ciotti 1991) is 


4^) = 




1 exp(—a;/t) 
_ 1 


dt, 


(B2) 


where s = r/Re and x = We have replaced the 

dummy variable t in Ciotti’s expression by 1/t, so as to 
avoid an integration to infinity. The physical value for the 
luminosity density is given by i/(r) = led’/Re u{s). Follow¬ 
ing Ciotti’s Figure 2, Figure Bl shows the luminosity density 
as a function of the normalised radius s, for a range of shape 
parameters n. 

The integrated spatial luminosity is then 


Ln(s)=47r f v„{s')s'^ds', 

Jo 


(B3) 


and is shown in Figure B2, with n ranging from 1 to 15. Here 


we have that L„(r) = led’R'e Ln{s), where Ln{r) yields the 
same total luminosity as given by equation p. We include 
Figures Bl and B2 in our paper for ease of reference and 
as they cover a greater range of shape parameters than in 
Ciotti’s Figure 2 (n = 2 — 10), where we encompass the 
n = 1 model being representative of an exponential disk 
and all integer values up to n = 15 being the cut oft used by 
Caon et al. (1993). 

For a spherical galaxy, with isotropic dynamics, con¬ 
stant M/L, and in hydrostatic equilibrium, the Jeans hy- 
drodynamical equation can be written as d{va^)/dr = 
G(M/L)L{r)v/d. Integration with the boundary condition 
va^ ^ 0 as r —> 00 gives 


\{s) = 


1 


el/s 


x(s) 


Ln{l/t)UnG/t)dt, 


(B4) 


where we have dropped the constants G{M/L). Again, we 
have replaced the dummy variable t with 1/t. This expres¬ 
sion for the velocity dispersion is dimensionless, the physical 
value u^{r) being given by 


ol{r) = G{M/L)Ied’Reol{s). 


(B5) 


The line of sight velocity dispersion is given by Ciotti 

1/77 


(1991) as 


niv) = 


LnG/t)Vn{l/t)yJl/t'^ — 772 ^ 4 ,(B 6 ) 


where 77 = R/Re- The luminosity weighted aperture velocity 
dispersion is then 


2 / N lo""" ^(V)^is{v)vdv 

<^ap,n(»7ap) - j^rj)ridr, ’ 


(B7) 


Figure Bl. The spatial luminosity density (ly) profiles are shown 
for a range in n from the exponential model, where n = 1, to 
n = 15, with s = r/Re being the 3D spatial radius divided by 
the projected half-light radius. (Extension of Fig. 2 from Ciotti 
(1991) which plotted values of n ranging from 2 to 10). 

Figure B2. The integrated spatial luminosity, normalised by the 
total profile luminosity, is plotted as a function of s = r//?e. 
(Extension of Fig. 2 from Ciotti (1991) which plotted values of n 
ranging from 2 to 10). 


where I{ri) is the profile. The physical quantity 

^ap,n(^ap) being given by the same transformation, snch 
that 

o'ap,n(-Rap) = G{M/L)Re’Re (ylp^„{rje,p)■ (B8) 

Shown in Figure are the spatial, line of sight, and aper¬ 
ture velocity dispersion as a function of radius for different 
profiles. It is noted that this quantity is dimensionless 
and should be multiplied by G{M/L)Iee^Re to obtain the 
physical quantity. The ratio of the spatial velocity dispersion 
at the spatial half-mass radius, s ~ 1.34, to the aperture ve¬ 
locity dispersion at the value of 77 corresponding to the ratio 
of the aperture size used and and the projected half-light 
radius of the galaxy, can be used to correct the measured 
central velocity dispersion to the physical value at the spa¬ 
tial half-mass radius. Thus, the multiplicative factor doesn’t 
come into play in this conversion process. 
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Table 1. The galaxy sample is from Bower, Lucey Sz Ellis (1992). Mulitple observa¬ 
tions (V-band) allow one to gauge the accuracy of the fitted model parameters due to 
various measurement errors and seeing conditions. Parameters from the best fitting 
and laws are shown. -Re,n is the effective half-light radius and Se,n the 

mean surface brightness within this radius. 


Sersic law de Vaucouleurs law (n=4) 


Galaxy 

Obs. 

n 

log i?e,n 

Se,n 

log Re,4 

Se,4 

Went. 

no. 


[kpc] 

[mag arcsec“^] 

[kpc] 

[mag arcsec 

N4339 

1 

4.26 

0.29 

20.89 

0.24 

20.73 


2 

4.74 

0.35 

21.13 

0.22 

20.65 

N4365 

1 

3.82 

0.43 

20.13 

0.47 

20.26 

N4371 

1 

3.24 

0.25 

20.11 

0.40 

20.64 


2 

3.10 

0.20 

19.96 

0.39 

20.62 

N4374 

1 

3.58 

0.34 

19.27 

0.43 

19.58 

N4377 

1 

5.35 

0.12 

20.10 

-0.03 

19.55 


2 

8.30 

0.46 

21.43 

-0.05 

19.46 

N4382 

1 

9.87 

1.79 

24.66 

0.54 

20.10 


2 

9.64 

1.62 

24.08 

0.42 

19.74 

N4387 

1 

2.36 

-0.19 

19.34 

-0.01 

19.98 


2 

2.17 

-0.20 

19.31 

-0.03 

19.90 

N4406 

1 

5.27 

0.89 

21.41 

0.58 

20.33 

N4435 

1 

4.43 

0.22 

19.75 

0.17 

19.56 


2 

3.80 

0.07 

19.22 

0.10 

19.32 

N4442 

1 

2.77 

-0.02 

18.68 

0.15 

19.28 


2 

2.96 

0.00 

18.73 

0.16 

19.33 

N4458 

1 

5.27 

0.41 

21.71 

0.21 

20.98 


2 

2.54 

0.01 

20.33 

0.19 

20.95 


3 

4.09 

0.18 

20.88 

0.16 

20.83 

N4464 

1 

2.33 

-0.40 

18.80 

-0.38 

18.80 


2 

2.63 

-0.41 

18.72 

-0.36 

18.89 

N4468 

1 

2.98 

0.08 

21.36 

0.27 

22.02 

N4472 

1 

3.31 

0.52 

19.58 

0.70 

20.19 

N4473 

1 

4.05 

0.09 

18.95 

0.08 

18.92 

N4476 

1 

3.28 

-0.10 

19.89 

-0.02 

20.18 


2 

3.06 

-0.13 

19.81 

-0.03 

20.19 

N4478 

1 

1.64 

-0.18 

18.70 

0.13 

19.80 

N4486 

1 

1.57 

0.20 

18.71 

1.00 

21.16 


2 

1.54 

0.19 

18.69 

1.01 

21.17 


3 

1.78 

0.25 

18.87 

0.94 

21.00 


4 

1.57 

0.21 

18.73 

0.97 

21.08 


5 

1.56 

0.21 

18.75 

0.99 

21.12 


6 

1.55 

0.19 

18.69 

1.07 

21.30 

N4486 

7 

1.66 

0.23 

18.78 

1.00 

21.15 


8 

1.61 

0.21 

18.74 

1.00 

21.16 

N4550 

1 

1.94 

-0.14 

19.07 

0.11 

19.94 


2 

1.85 

-0.16 

19.02 

0.17 

20.19 

N4551 

1 

2.24 

-0.10 

19.61 

0.16 

20.51 

N4552 

1 

3.99 

0.18 

19.03 

0.19 

19.04 


2 

4.12 

0.22 

19.17 

0.20 

19.10 


3 

4.30 

0.24 

19.25 

0.19 

19.07 


4 

5.44 

0.41 

19.87 

0.19 

19.06 


5 

5.12 

0.36 

19.69 

0.19 

19.06 


6 

6.46 

0.56 

20.43 

0.20 

19.08 


7 

5.36 

0.42 

19.92 

0.23 

19.21 


8 

5.45 

0.40 

19.85 

0.19 

19.07 


9 

4.44 

0.27 

19.34 

0.20 

19.08 


10 

4.96 

0.32 

19.53 

0.18 

19.01 


11 

4.12 

0.21 

19.13 

0.19 

19.06 


12 

5.18 

0.37 

19.74 

0.20 

19.11 


13 

3.96 

0.17 

18.99 

0.18 

19.02 


14 

5.70 

0.43 

19.95 

0.18 

19.05 

N4564 

1 

3.82 

0.09 

19.59 

0.11 

19.68 


2 

5.47 

0.25 

20.21 

0.07 

19.51 

N4621 

1 

4.01 

0.43 

20.00 

0.43 

20.00 


2 

5.90 

0.65 

20.83 

0.30 

19.55 
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Table 1. - continued 

Sersic law de Vaucouleurs law (n=4) 


Galaxy 

Obs. 

n 

log Re,n 

Se,n 

log Re,4 

Se,4 

Ident. 

no. 


[kpc] 

[mag arcsec“^] 

[kpc] 

[mag arcsec 

N4636 

1 

3.63 

0.57 

20.61 

0.66 

20.93 


2 

3.78 

0.64 

20.86 

0.70 

21.04 

N4660 

1 

3.41 

-0.19 

18.46 

-0.15 

18.60 


2 

3.49 

-0.19 

18.44 

-0.15 

18.58 


3 

3.67 

-0.19 

18.44 

-0.16 

18.56 


4 

3.54 

-0.20 

18.41 

-0.17 

18.53 

N4697 

1 

4.64 

0.72 

20.76 

0.58 

20.26 
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Table 2. The galaxy sample is from Bower, Lucey Sz Ellis (1992). Parameters from the 
best fitting and profiles have been corrected to a form suitable for use in the 

construction of the Fundamental Plane. Re ,4 and Re,n are the spatial (3 dimensional) 
half-light radii of the R^^^ law and the law respectively, with Se ,4 and Se,n the 

mean surface brightness term within R^. ctq is the central aperture velocity dispersion 
measure, crj^/ 2,4 ^ 1 / 2 ,n spatial velocity dispersions at the spatial half-light 

radius of the and the R^^^ profiles respectively (as detailed in the text), and crtot,n 

is the aperture velocity dispersion, for the model, with radius extended to infinity. 


Galaxy 

o-Q 

n 

logi?e,n 

Se,n 

<^ 1 / 2 ,n 

crtot,n 

log Re ,4 

Se,4 

0'l/2,4 

Ident. 

km 


[kpc] 

[mag] 

km 

km 

[kpc] 

[mag] 

km 

N4339 

112 * 

4.50 

0.45 

21.00 

79 

84 

0.36 

20.69 

83 

N4365 

258 

3.82 

0.56 

20.13 

188 

195 

0.60 

20.26 

185 

N4371 

127 

3.17 

0.35 

20.03 

100 

101 

0.53 

20.63 

92 

N4374 

302 

3.58 

0.47 

19.27 

227 

233 

0.56 

19.58 

217 

N4377 

136 

6.82 

0.46 

20.57 

83 

97 

0.09 

19.50 

107 

N4382 

200 

9.76 

1.85 

24.33 

72 

94 

0.61 

19.91 

143 

N4387 

115 

2.26 

-0.07 

19.32 

103 

102 

0.11 

19.94 

90 

N4406 

226 

5.27 

1.02 

21.41 

138 

152 

0.71 

20.33 

160 

N4435 

171 

4.11 

0.28 

19.45 

128 

134 

0.26 

19.43 

130 

N4442 

217 

2.87 

0.12 

18.70 

181 

182 

0.29 

19.30 

163 

N4458 

89 

3.97 

0.36 

20.84 

67 

70 

0.32 

20.92 

67 

N4464 

120 

2.48 

-0.28 

18.76 

112 

111 

0.24 

18.84 

102 

N4468 

74* 

2.98 

0.21 

21.36 

61 

61 

0.40 

22.02 

55 

N4472 

298 

3.31 

0.65 

19.58 

225 

229 

0.83 

20.19 

209 

N4473 

185 

4.05 

0.22 

18.95 

141 

148 

0.21 

18.92 

142 

N4476 

79* 

3.17 

0.02 

19.85 

67 

68 

0.11 

20.18 

62 

N4478 

148 

1.64 

-0.05 

18.70 

136 

133 

0.26 

19.80 

112 

N4486 

337 

1.61 

0.34 

18.74 

294 

287 

1.13 

21.14 

236 

N4550 

81* 

1.89 

- 0.02 

19.04 

74 

72 

0.27 

20.06 

61 

N4551 

105 

2.24 

0.04 

19.61 

92 

91 

0.29 

20.51 

79 

N4552 

246 

4.90 

0.47 

19.49 

169 

182 

0.32 

19.07 

184 

N4564 

153 

4.64 

0.31 

19.86 

110 

118 

0.22 

19.59 

117 

N4621 

218 

4.95 

0.68 

20.34 

143 

155 

0.50 

19.75 

158 

N4636 

201 

3.70 

0.74 

20.73 

146 

150 

0.81 

20.98 

141 

N4660 

183 

3.53 

-0.06 

18.44 

153 

157 

0.03 

18.57 

148 

N4697 

189 

4.64 

0.85 

20.76 

124 

133 

0.71 

20.26 

133 


* Corrected value taken from McElroy (1995). 
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Table 3. Fractional Variance from the PCA Analysis 


model 

Vari 

Var2 

Vara 

O-Q 

56.47% 

42.26% 

1.27% 


63.68% 

35.68% 

0.64% 

0-1/2,n 

64.78% 

34.40% 

0.82% 


61.96% 

37.30% 

0.74% 
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Table 4. Fundamental Plane R oc via the 3D bisector 

method 


model 

R oc 


rms 


A 

B 

residual 


1.33±0.10 

-0.79±0.11 

0.053 

-Rl/",cro 

1.26±0.10 

-0.78±0.08 

0.053 


1.48±0.13 

-1.00±0.09 

0.062 


1.44±0.11 

-0.93±0.08 

0.057 
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Table 5. Regression Analysis for log Re versus (logcr + 6Ee). All results are from the 
program SLOPES (Feigelson Babu 1992). The reader is referred to their paper for 
an understanding of each method of analysis. The approach referred to as uses 

the projected half-light radius from the model, the mean surface-brightness within 

this radius and the projected central velocity dispersion in the construction of the FP. 
R^/”'[(To], R^^”'[(Ttot,n] use the spatial half-light radius from the 

model and the associated surface brightness term which allows for different luminosity 
structures in each galaxy as represented by n. These approaches use the projected central 
velocity dispersion term, the spatial half-light velocity dispersion term and the infinite 
aperture velocity dispersion term respectively. 


method 

Asymptotic Formulae 

Bootstrap 

Jackknife 


Slope (A) 

Slope (A) 

Slope (A) 

_Rl/4[o-(,] . f, = 0.236 ± 0.031 

OLS(Y|X) 

1.256±0.080 

1.257±0.088 

1.256±0.088 

OLS(X Y) 

1.377±0.067 

1.384±0.076 

1.378±0.079 

OLS bisector 

1.315±0.069 

1.318±0.077 

1.315±0.077 

Orthogonal 

1.332±0.070 

1.336±0.080 

1.332±0.079 

Reduced major axis 

1.315±0.069 

1.319±0.077 

1.315±0.077 

Mean OLS 

1.317±0.068 

1.321±0.076 

1.317±0.077 

b = 0.254 ±0.027 

OLS(Y|X) 

1.200±0.050 

1.198±0.064 

1.200±0.058 

OLS(X Y) 

1.249±0.044 

1.256±0.057 

1.250±0.052 

OLS bisector 

1.224±0.045 

1.226±0.058 

1.224±0.049 

Orthogonal 

1.229±0.045 

1.232±0.059 

1.230±0.050 

Reduced major axis 

1.224±0.045 

1.226±0.058 

1.225±0.049 

Mean OLS 

1.225±0.045 

1.227±0.058 

1.225±0.049 

b = 0.285 ± 0.032 

OLS(Y|X) 

1.347±0.073 

1.332±0.090 

1.346±0.100 

OLS(X Y) 

1.439±0.054 

1.442±0.071 

1.439±0.058 

OLS bisector 

1.392±0.059 

1.385±0.074 

1.391±0.072 

Orthogonal 

1.407±0.058 

1.403±0.075 

1.407±0.067 

Reduced major axis 

1.392±0.059 

1.386±0.074 

1.392±0.072 

Mean OLS 

1.393±0.059 

1.387±0.073 

1.392±0.071 

: b = 0.270 ± 0.029 

OLS(Y|X) 

1.333±0.067 

1.324±0.083 

1.333±0.086 

OLS(X Y) 

1.408±0.053 

1.413±0.068 

1.409±0.058 

OLS bisector 

1.370±0.057 

1.367±0.071 

1.370±0.065 

Orthogonal 

1.382±0.056 

1.381±0.072 

1.382±0.062 

Reduced major axis 

1.370±0.057 

1.368±0.071 

1.370±0.065 

Mean OLS 

1.371±0.056 

1.369±0.071 

1.371±0.064 
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